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ABSTRACT 


The realization of linear-phase-shift and constant-time-delay trans- 
fer characteristics using computer-aided analysis and synthesis of a 
fourth-order highpass active filter was the main objective of this investi- 
gation. 

Two approaches were to be taken in the analysis and synthesis phase 
Only the one using topologically oriented programs was completed. The 
second involved the derivation and analysis utilizing the system trans- 


fer function. 
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I. INTRODUCTION 


The primary objective of this analysis was to determine the parameter 
modifications to a Sallen-Key type RC active filter to give linear-phase- 
shift and constant-time-delay transfer characteristics. In 1955 Sallen 
and Key [Ref. 1] published a catalog of RC active filters. The circuit, 
Fig. 1, chosen for investigation is a solid-state version of one such 
filter. The chosen circuit was designed by O. M. Baycura, Associate 
Professor of Electrical Engineering at the Naval Postgraduate School. 

The network is actually a fourth-order highpass active filter made up by 
cascading two second-order filter sections. 

Two approaches to the analysis were to be taken, both of which 
solicit the aid of a digital computer. One method involved network analy- 
sis using a topology-oriented program such as ECAP. The second involved 
the derivation and investigation of the system transfer function. Only 
the first approach to the problem was completed, due to difficulties en- 
countered with the computer program, and will be the primary topic of 
discussion. 2 

Upon completion of the computer-aided analysis phase, the filter 
was constructed and tested using the determined parameters. The actual 
filter characteristics and those calculated by the computer program were 


then compared. 
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II. CIRCUIT DESIGN 


As previously discussed, the circuit is a fourth-order highpass RC 
active filter made up of two cascaded second-order sections. The second- 
order filter network is considered by some as the standard building block 
in filter design. The obvious advantage of : second-order system over 
higher-order eyerens is in the simplicity of analysis and optimization. 
Higher-order filter networks can be realized by cascading individually 
tuned sections. Figure 2 shows the basic building block of the active 
filter that-was investigated. This circuit configuration is described by 
Mullaney [Ref. 2] as an RC peaking network and is one of the most 
widely used filter types. 

Assuming the active devices are unity-gain amplifiers, the network 


is described in Ref. 2 by the transfer function: 





K(s/uy,)? 
a (1) 
(s/w) + s/Q wit 
where, 
W = response break~frequency in radians 
oe: 
Q= on ° quality factor = amount of peaking of 


, a the bandpass filter. 
The case with unity-gain amplifiers is a special one and the general case 
is one in which the amplifier gain is represented by "A". The transfer 
function describing this case is derived using Fig. 2 with the gains 
of amplifiers given by A, and A, . These amplifiers are assumed to be 


ideal and represent dependent sources whose input and output 
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General Second-Order Highpass Filter 
Figure 3. 
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characteristics do not affect the circuit. This leaves the network with 
the passive complement of two resistors and two capacitors which will 
result in a second-order polynomial in the denominator of the transfer 

function. The general form of the transfer function was derived in the 

following manner using the circuit in Fig. 3: 


Starting at Node 5, 


ent : Dy a. V3 (2) 
where, 
: S 
i ee )*V, (3) 
272 
Vo = A, x Vy (4) 
a = ee. 
ol 7 Vin = ( s + VA Ga ) id Voue: (5) 


Substituting equations (3) and (4) in (2) and solving for Me in terms of 


V 
out 
s+1/R.C V 
(ee Houta \ = 
My et S a a ALA, ) eS (6) 


Substituting (6) in (5) and solving for the ratio. of The to Vi! 


ee 


V 
; Yd , 
in S +sl1/R,C,(1 -A,)+1/R,C, J+ 1/2, CRC, 

Let 


W, = 1/R,C and W, = 1/R,C 


] Z Z 


Wo = fw, w, . 
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Then equation (7) becomes 


A Ba 
G(s) = | —— 2? (9) 
[ s+s[w (1-A)+w,] + we | 


where the quality factor is 


JR,C,R,C, 
RoC, - A) +R Cy 


The transfer function derived reduces to the form in equation (1) when 


the amplifier gains, A, and Ay are set to unity. For convenience, the 


1 
denominator of the transfer function was normalized in the following 


manner. 
D(s) = (s/w) + (s/w )d + 1 


where 
d= 1/@ = the dampime factor. 
Analysis of second-order systems shows that this configuration is 
characterized by its two conjugate poles anda pair of zeros at the ori- 


gin. The roots of the denominator can be found by using the quadratic 


JHda/w Ne - ave 
O O 


2 


formula. The roots are given by 
s=-d/2w + 
’ O = 


which shows that they lieon a circle of radius Wo and have real parts 
equal to -d/2) |. The system is unstable, (the poles lie in the right- 


hand side of the s-plane) when the discriminant is positive real and 
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The nominal value for the measured voltage gain of the emitter- 
follower amplifier used was 0.6 at 5 Mhz, and assuming the gains A, and 


A, equal, a value for the damping factor of 1.7 is obtained from the graph. 
Using the hybrid-pi model of the transistor, and the analysis procedures 
discussed in the text by Millman and Halkias [Ref. 47) , an attempt was 
made to calculate the theoretical voltage gain. Difficulties were en- 
countered at this point in determining the necessary parameters. First, 
published data for the transistor type used was unavailable and, secondly, 
attempts to measure the high-frequency hybrid parameters failed to pro- 
duce consistent and useful results. The decision was made to continue 
the investigation by specifying the parameters for a suitable active device. 


The hybrid-pi model as discussed in Ref. 4 is shown in Fig. 5 in its 


high-frequency configuration. 


Chic 
B Tb 'b pte c 
e 6 
Thte Chte Em’ bte “ce 
E 
e @ @ 


Figure 5. Hybrid-pi Equivalent for Transistor in CE Configuration 


After some research on the subject, a model of a transistor studied 
by Leon [Ref 5. ] was chosen as a guideline for specifying the assumed 


parameters. These parameters are listed in Table 1. With the equivalent 
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circuit for the active device selected, the equivalent circuit in Fig2 6 


for the fourth-order active filter was completed. This is the circuit as 


initially programmed in the topological approach of the analysis. 


Nagnitude 


Base-spreading Resistance r 

Input Resistance r 

Output Resistance r 

Reverse Feedback Resistance r 
Emitter-base Junction Capacitance C 


Collector Junction Capacitance C 


Transistor Transconductance g 





60 ohms 

lel) kilohms 
25 kilohms 
10 megohms 
38 picofarads 
3 picofarads 


90 millimhos 


Table 1. Hybrid-pi Parameter Values 
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Ill, COMPUTER-AIDED ANALYSIS 

At this stage, an existing computer program or programs to aid in 
the analysis of the network was selected. The topological approach to 
the investigation required that the program be capable of accepting topo- 
logical circuit description as an input. Since a lumped model with the 
parameters independent of frequency was used as the active device's 
equivalent, the program need only be able to solve linear or small-signal 
transistor circuits. The output capabilities necessary were the following: 

Calculate the frequency response and make a Bode plot. 

Determine the time delay as a function of frequency and plot. 

Solve for the transfer function of the circuit. 

Calculate and plot the transient response of the system. 

Being restricted to the program library maintained by the W. C. Church 
Feompiter Center at the Naval Postgraduate School, the selection of a 
suitable program was narrowed considerably. 

An integrated package of Fortran IV programs called LISA, for Linear 
Systems Analysis, was located. LISA fulfilled all of the requirements 
except for the plot of the time delay. The inputs to LISA as described 
in the user's guide [Ref. 6] may be a topological circuit description, 
matrix equation, transfer function, or block diagram. The program uses 
Laplace transform techniques for analysis and the user is allowed to 
work in the s-plane, or in the frequency and time domain. The topology 
and the component values of the circuit can be altered by simple pro- 


gramming and they are outputed for comparison and examination. 


uS 





To use LISA in the analysis of the active filter, the equivalent cir- 
cuit description is entered in a preformated form as values of R, L, C, 
mutual inductance, initial conditions, independent current sources, and 
voltage- or current-dependent current sources. The dependent sources 
are utilized in modeling active devices in the hybrid-pi or Tee equiva- 
lent transistor circuits. The program then forms the nodal matrix equation 
in the complex frequency variable, s=o+jwW. Also, the functions re- 
quired by the problem must be defined. These can be voltages, voltage 
ratios, currents, current ratios, admittances, impedances, determinants, 
co-factors, or any three-terminal, grounded, two-port network parameters. 
The program can compute and output as requested poles and Zeros, a root 
locus, frequency and transient response, or sensitivity. The program is 
capable of computing the Bode response and sensitivity of circuits of 
up to a maXimum of 50 nodes. LISA is also programmed to output error 
messages for errors detected in the input and compute phases. 

After the nodal matrix has been formed, the program is ready to 
enter the COMPUTE phase. There are three wemrees by which LISA solves 
the matrix equations. 

The first option, named POLY, solves for the polynomials that define 
the poles and Zeros of the requested functions. POLY can handle systems 
or order 12 or less with proper input scaling. The limitation of this 
option is the round-off errors and overflow during floating-point calcu- 
lations that is caused by its inability to cope with magnitude differences 
greater than 10° ; 
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The second method is TRFN which finds the roots of the function 
directly using the determinant that defines them. The limitations found 
in POLY are avoided by TRFN since programmed scaling is done at critical 
locations in the algorithm. The maximum size system sucessfully solved 
DyeuRE IN is; 33. 

The third option is ACCA which works in the real frequency domain 
only. ACCA determines the system response for s = jw by evaluating the 
polynomials in the matrix equation and obtains a system of linear equa- 
tions with complex coefficient matrix. The system of equations is solved 
using a modified method of Gaussian elimination. The size of the sys- 
tem that can be handled by ACCA is limited only by the storage available. 
Systems of the order 50 have been solved by ACCA. 

Because of these three options it is possible to check the accuracy 
of one against another by solving the same problem with each. The out- 
puts of each can be in the form of a listing or can be plotted. The results 
of the POLY and TRFN options are polynomials and/or roots. To geta 
frequency-response listing or plot, the option BODE must be used. The 
ACCA method of solution results in the computation of magnitude and 
phase over a given range of frequency and the output can be a listing 
and/or plot. 

In addition to the magnitude and phase being computed, the system 
time delay is also computed. The time delay is found by differentiating 
the phase function with respect to frequency; that is, 


Time delay t(w) = oF 6 (w) ° 


Z\ 





The output of the time delay is not plotted, but is listed by all three 
when the Bode plot is requested. 

The equivalent circuit in Fig. 6 was programmed on LISA without any 
scaling. The network consists of one ground node and 16 circuit nodes. 
The hybrid-pi parameters were renamed for convenience and an independent 
current source and source resistance were added. The TRFN option for 
solving the system was chosen and the BODE option, used to obtain the 
frequency response in the range from 10 Hz to 10 GHz, was requested 
as a listing and a Bode diagram. The impulse response of the system was 
requested for the first 150 nanoseconds. The listing and the plots gen- 
erated could not be included as received due to their large size. The 
listings will be omitted and the plots, which are made on the line printer, 
redrawn for proper size for inclusion in this report. 

The method of attack followed in the analysis was that of the 
cascading of individually tuned sections to achieve the desired results. 
The first second-order filter section was left untouched while the second 
Pa 


section was tuned. This meant changing the values of R ; Ros and 


] ] 


C,. which correspond to R respectively, Fig. 6. 


9 C,, R,, and C 


3’ ~3 


The results of the computation with network topology as specified 


& 


originally were disappointing. Tuning the second section did not pro- 
duce the results desired, namely linear phase shift and constant time 
delay. The fourth-order system was not actually of order four due to the 
junction capacitances assumed in the active device model. In the analy- 
sis done by Sallen and Key discussed in Ref. 1, the active devices are 


Z2 





assumed to be ideal. The amplifiers are assumed to have infinite input 
impedance, Zero output impedance and stable gain. These conditions are 
also assumed by Mullaney [Ref. 2]. Optimization with the circuit con- 
figured as in Fig. 6 seemed futile if not unrealistic at this point. Thus, 
a decision was made to simplify the problem by eliminating the junction 
capacitances, programming the network as revised in Fig. 7. 

The voltage gain of the amplifier can be Paine using the pro- 
cedures and relations developed in Ref. 4. Converting the hybrid-pi 
parameters in Table 1 to hybrid parameters, the voltage gain computes to 
0.945 giving a damping factor from graph of Fig. 4 of 0.84. The list of 
the input statements is given in Table 2 el the frequency and time re- 
sponses in Fig. 8. In Fig. 8, the values of the resistors R, and R, are 


1 2 
3.3 kilohms and 6.6 kilohms respectively and the capacitors C, and Cy 
are equal at 10 picofarads. 


The first phase of the search for the optimum values of capacitors, 


C,-ane’C 


3 4: was carried out using the ACCA method since it possessed 


its own scaling feature. Note that the magnitudes of the components 
programmed go from a high of 6.6 kilohms to a low of oie farads. 

Even with its own scaling program, ACCA encountered difficulties in com- 
putation. Negative time delays were found computed on several runs. 
There was no correlation of this problem with the order-of-magnitude 
difference between Cy and Cy. and C, and Cy: Using ACCA, negative 
time delays were computed when the values of all the capacitors were set 


equal to 10 pf and then one pair, C, and Cys was held at 10 pf and the 


23 


= i =. 





qno 


juaTeainby Aouenbelj MOT - Je UII satpy ssedyhHty 1epiO-yynoj 


°74 oiInbii 


0 





24 





Table 2 LISA Input Statements 


TITLE, FORTH-ORDER HIGHPASS ACTIVE FILTER 
NOTE, HYBRID-PI TRANSISTOR EQUIVALENT CIRCUIT 
READ, CIRCUIT 


RB1, 2.360 

REL, 3,4, 1.483 

R34, 3,500 

TGM1, 5,351, O9E=3 

RL1, 5,0, 100 

C2,4,6, 1E-11, TOL=.05 

R2,6,0,6.653 

RB2,6,7,60 

RD2,7,8,1.4E3 

Ri, 8,0,500 

TGM2, QyTy 8, eO9E=3 

RL2, 9,0, 100 

Ra, 8, 2.34383 

oe: ic alo pee ala 

RB3,10,11, 60 

MOS vil, be, .uh3 

R7,12,0, 500 

TGM3,13,11,12, .09E-3 

RL3, 13,0, 100 

Cri) Qe, = 1 

R6,14,0,6.6E3 

RBY, 14,15, 60 

RD, 15,16, 1.453 

R8,16,0,500 

TGM1,17,15,16, .09E-3 

Riven «©. 100 

R5,16,0, 3.3E3 

DEFINE, VOUT=V(16) 

COMPUTE, TREN, VOUT 

DATA, FREQ=20, 10, 10.0, Hz 

COMPUTE, BODE, VOUT 

PPLOT 

LABEL, C1=C2=C3=Ch=10 pf 

READ, DRIVER 

IMPULSE, 0,10. 

DATA, TIME=0,1.5E-7, 1E-9 

COMPUTE, TRANSIENT, VOUT 

PPLOT, TIME, VOUT 

LABEL, VOUT FOR IMPULSE 

EXIT 

Frequency Range: 10 Hz to 10 GHz in 10 decades with 20 - 
points per decade. 


Time Interval: 0 to 150 nanoseconds at 1 point per ns. 
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other pair changed to 2 pf. When this phenonemon occurred, the nega- 
tive times were confined to the region well below the passband. Also, 
associated with the negative times were output magnitudes of less than 
10°”. 

The TRFN method of solution in combination with the BODE option was 
next checked for this error. Negative time delays were found to occur 
when using this method also, but not in all the cases that were encountered 
using ACCA. Still no correlation could be found between capacitor magni- 
tudes and the phenomenon. Since TRFN has no programmed scaling fea- 
ture, several runs were repeated with the component magnitudes scaled. 
Improvement was noticed, but it was not overwhelming. The negative time 
delays were found to occur with about the same pattern seen in the ACCA 
runs. 

The conclusion was that there must be some order-of-magnitude 
limit that the program can handle successfully. Just as with the option 
POLY, the subroutines could be bothered by the same limitations that 
plague the former. Many root-finding programs have difficulties when 
the order of magnitudes of the roots becomes very large or the time con- 
stants in the system differ by many orders of magnitudes. This is what 
most likely occurred in this analysis. 

Eliminating those parameter values that are troublesome, the range 
of optimumization was narrowed to between one and two picofarads. The 
value at the iene end, 1 pf, results in computation of negative time 
delays. No difficulties were encountered using any higher values. The 
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results of the last phase topological analysis approach are shown in 


Figs. 9 and 10 with the values determined listed in Table 3. 
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Figure 9. Impulse Response 


Filter Time Delay 
Magnitude Cutoff Passband 
C3 Ch, Frequency f, Attenuation maximum 
pf MHz : dbs ns 


2.73x102 2.7x107° 


2.6lyx102 3elx107! 
9.56x102 heox1o~! 


negative h.7x107/ 


Table 3. Results of Computer-Aided Analysis - Topological Approach 
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Figure 10. Frequency Response with Optimum Parameter and Actual Circuit (@) 





IV. CIRCUIT PERFORMANCE 


The original circuit, Fig. 1, less the output amplifier, was modified 
by changing the values of Cy and C, to 2 picofarads. This value is the 
upper bound on the range of optimum values determined in the analysis. 

The frequency Peseeeriiae of the filter was then taken. The results 
were quite different than predicted and are plotted in Figure 10. The 
attenuation through the actual filter was on the average 10 to 15 db greater 
outside the passband. There was a definite downward slope in the pass- 
band. The cutoff frequency was also shifted to almost half of its calcu- 
lated value. 

A first the measurement techniques and instrumentation were suspect, 
but this proved to be false after careful review of procedures and equip- 
ment operation. Since the circuit was constructed ona "bread-board", 
stray and wiring capacitances were next thought to be the causes, and 
the circuit was transferred to printed circuit form to eliminate them. The 
frequency response was taken again. The results did not change enough 


to warrant another plot. 


and C,, were tried and 


Next, other values for the capacitors, Co 4 


the only significant changes noticed were outside the passband. 

The minor adjustments suggested by Sallen and Key in Ref. 1 did 
not help to achieve the designed performance. The gains of the ampli- 
fiers were next altered to see their effect on the response characteristic. 
The gains were increased slightly by removing the collector resistors 
and decreased by increasing the collector resistance. There was some 
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improvement in the shape of the characteristics when the gain was 
increased. 
The problem, it appears, goes a little deeper than the physical com- 


ponents and their minor adjustments. 
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V. CONCLUSIONS 


The circuit would have performed as predicted if the active network 
constructed with physical components has possessed the assumed ideal 
parameters: infinite input impedance, zero output impedance, and highly 
stable gain. Referring to equations (9) and (10) on page 14, the damping 
factor d is directly related to the gains of the active devices. The damp- 
ing factor influences the shape of the frequency response characteristics. 
The next step, if the analysis were to continue, would be to replace 
the junction transistors with field-effect types. The field-effect tran- 
sistor is noted for its high input impedance. The output impedance is no 
greater than that of a junction transistor. The capacitances associated 
with the FET's are also small and no greater than those found in junction 
transistors. 

An error in judgement was committed when the decision was mde to 
continue the analysis and synthesis at high frequencies without consider:- 
ing the effects of the junction capacitances. Had these factors been 


taken into account, the synthesis could have been more successful. 
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